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Abstract 

^ ' We introduce a spectral approach to non-perturbative field theory 

l-l , within the periodic field formalism. As an example we calculate the real 

and imaginary parts of the propagator in 1 + 1 dimensional (^* theory, 
identifying both one-particle and multi-particle contributions. We dis- 
cuss the computational limits of existing diagonalization algorithms and 
r^ ' suggest new quasi-sparse eigenvector methods to handle very large Fock 

O^ , spaces and higher dimensional field theories. [PACS numbers: ll.lOKk, 

12.38Lg] 

^. 1 Introduction 

^ ■ 

Modal expansion methods have recently been used to study non-perturbative 
phenomena in quantum field theory ||l||-P|. Modal field theory, the name for 
K^ , the general procedure, consists of two main parts. The first is to approximate 

\^ ' field theory as a finite-dimensional quantum mechanical system. The second 

is to analyze the properties of the reduced system using one of several com- 
putational techniques. The quantum mechanical approximation is generated 
by decomposing field configurations into free wave modes. This technique has 
been explored using both spherical partial waves (spherical field theory ||l|l@) 
and periodic box modes (periodic field theory Q). 

Having reduced field theory to a more tractable quantum mechanical sys- 
tem, we have several different ways to proceed. Boson interactions in Euclidean 
space, for example, can be modeled using the method of diffusion Monte Carlo. 
In many situations, however, Monte Carlo techniques are inadequate. These 
include unquenched fermion systems, processes in Minkowski space, and the 
phenomenology of multi-particle states. Difficulties arise when the functional 
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integral measure cannot be treated as a probability distribution or when in- 
formation must be extracted from excited states obscured by dominant lower 
lying states. Fortunately there are several alternative methods in the modal 
field formalism which avoid these problems. Matrix Runge-Kutta techniques 
were introduced in pj as a method for calculating unquenched fermion interac- 
tions. Here we discuss a different approach, one which directly calculates the 
spectrum and eigenstates of the Hamiltonian. For this approach it is essen- 
tial that the Hamiltonian is time-independent, and so we will consider periodic 
rather than spherical field theory. As we demonstrate, these methods naturally 
accommodate the study of multi-particle states and Minkowskian dynamics. 

We apply the spectral approach to H- 1 dimensional (/)^ theory in a periodic 
box and calculate the real and imaginary parts of the 4> propagator. Some in- 
teresting properties of (/)| theory such as the phase transition at large coupling 
were already discussed within the modal field formalism using Euclidean Monte 
Carlo techniques Q]. The purpose of this analysis is of a more general and 
exploratory nature. Our aim is to test the viability of modal diagonalization 
techniques for quantum field Hamiltonians. We would like to know whether 
we can clearly see multi-particle phenomena, the size of the errors and com- 
putational limitations with current computer resources, and how such methods 
might be extended to more complicated higher dimensional field theories. 

The spectral method presented in the first part of our analysis is similar 
to the work of Brooks and Frautschi H],R who considered a 1 + 1 dimensional 
Yukawa model in a periodic box and deserves credit for the first application of 
diagonalization techniques using plane wave modes in a periodic box. Our cal- 
culations are also similar in spirit to diagonalization-based Hamiltonian lattice 
formulations B and Tamm-Dancoff light-cone and discrete light-cone quanti- 
zation [p|-[pl. There are, however, some differences which we should mention. 
As in ff we are using a momentum lattice rather than a spatial lattice. We 
find this convenient to separate out invariant subspaces according to total mo- 
mentum quantum numbers. Since we are using an equal time formulation our 
eigenvectors are not boost invariant as they would be on the light cone. Also we 
are using a simple momentum cutoff scheme rather than a regularization scheme 
which includes Tamm-Dancoff Fock-space truncation. As a result our renor- 
malization procedure is relatively straightforward, but we will have to confront 
the problem of diagonalizing large Fock spaces from the very beginning. In the 
latter part of the paper we mention current work on quasi-sparse eigenvector 
methods which can handle even exceptionally large Fock spaces. Despite the 
differences among the various diagonalization approaches to field theory, the 
issues and problems discussed in our analysis are of a general nature. We hope 
that the ideas presented here will be of use for the various different approaches. 
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2 Spectral method 

The field configuration <^ in 1 + 1 dimensions is subject to periodic boundary 
conditions 0(i, x— L) ~ (p{t, x + L). Expanding in terms of periodic box modes, 
we liave 



n=0,±l,... 

The sum over momentum modes is regulated by choosing some large positive 
number iVmax and throwing out all high- momentum modes 0„ such that \n\ > 
-^max- In this theory renormalization can be implemented by normal ordering 
the cf)'^ interaction term. After a straightforward calculation (details are given 
in 1^), we find that the countertcrm Hamiltonian has the form 



Tl — — _ 

where 



6A6 

4I2Z- 



v-^ 1 /.,2_2 



b= 1^ ik^ ^n = ^'^+P?. (3) 

n=-Ar„ax,Wmax 

We represent the canonical conjugate pair (/)„ and TJ" using the Schrodinger 
operators </„ and —i-^—. Then the functional integral for cj)'^ theory is equivalent 
to that for a quantum mechanical system with Hamiltonian 

H= E [-I0h:^. + '2i-l-m<l-n^n\ (4) 

n = -JV„»x,Wmax 

ni = -iV„ax,iVmax 
"1 +"2 +"3 +"■4=0 

We now consider the Hilbert space of our quantum mechanical system. 
Given d, an array of non- negative integers, 

d = {d-Ar„„, , ■ • ■ rfjv,„ax} ' (5) 

we denote Pdil) as the following monomial with total degree \d\, 

pM= n 9""' E'^-^l'^l- (6) 

n=-JVmax,JVmax " 

We define G(^{q) to be a Gaussian of the formg 



Gciq) = Yi ^^P 



n=-N,j,^^.N„ 



2 



(7) 
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G(^{q) has been defined such that Gfj.{q) is the ground state of the free theory. 



C is an adjustable parameter which we will set later. Any square-integrable 
function ip^q) can be written as a superposition 



i>iq)=J2cdPdiq)Gc{q). (8) 

d 

In this analysis we consider only the zero-momentum subspace. We impose this 
constraint by restricting the sum in (@) to monomials satisfying 



y^ndn 



0. (9) 



We will restrict the space of functions ipiq) further by removing high energy 
states in the following manner. Let 

fc(d)=^|n|d„. (10) 

n 

k{d) was first introduced in pi and provides an estimate of the kinetic energy as- 
sociated with a given state. Let us define two auxiliary cutoff parameters, -R'max 
and I?max- We restrict the sum in (m to monomials such that k{d) < K^aax 
and |d| < Umax- We will refer to the corresponding subspace as VK^^^,D^a^- 
The cutoff ii'max removes states with high kinetic energy and the cutoff -Dmax 
eliminates states with a large number of excited modes.R We should stress that 
-K^max and -Dmax are only auxiliary cutoffs. We increase these parameters until 
the physical results appear close to the asymptotic limit -ftTmax, -Dmax ^ oo. In 
our scheme ultraviolet regularization is provided only by the momentum cutoff 
parameter ./Vmax- 

Our plan is to analyze the spectrum and eigenstates of H restricted to this 
approximate low energy subspace, Vk^^^,d^^^- For any fixed L and ./Vmax, H is 
the Hamiltonian for a finite-dimensional quantum mechanical system and the 
results should converge in the limit -R'max, -Dmax —^ oo- We obtain the desired 
field theory result by then taking the limit L, ^^s^ -^ oo. 



3 Results 

We have calculated the matrix elements of H restricted to Vk^^^,d^^^ using a 
symbolic differentiation-integration algorithmn and diagonalized H, obtaining 
both eigenvalues and eigenstates. Let A be the full propagator, 

A(p2) =i ( d^x e^P-^" (0| T [(j){x^')(j){Q)] |0) . (11) 



•^In the case of spontaneous symmetry breaking, the broken symmetry of the vacuum may 
require retaining a large number of go modes. This however is remedied by shifting the 
variable, q'g = qq - (0| ijo |0>. 

■^AU codes can be obtained by request from the authors. 



We have computed A by inserting our complete set of eigenstates (complete in 
^ifmax.i'max)- ^ct A^p bc thc multi-particlc contribution to A, 

A,„p(p2) = A(p2)_Ap„i^(p2)^ (12) 

where Apoic is the single-particle pole contribution. We are primarily interested 
in Amp, a quantity that cannot be obtained for p^ > using Monte Carlo 
methods. Since the imaginary part of Apoio is a delta function, it is easy to 
distinguish the single-particle and multi-particle contributions in a plot of the 
imaginary part of A. The real part of A, however, is dominated by the one- 
particle pole. For this reason we have chosen to plot the real part of Amp rather 
than that of A. 

Although we have referred to multi-particle states, it should be noted that 
in our finite periodic box there are no true continuum multi-particle states. 
Instead we find densely packed discrete spectra with level separation of size ^ 
L^^ which become continuum states in the limit L —> oo. We can approximate 
the contribution of these L — )■ oo continuum states by a simple smoothing 
process. We included a small width F ~ L~^ to each of the would-be continuum 
states and averaged over a range of values for L. For the results we report here 
we have averaged over values L — 2.0tt, 2.1tt, ■ ■ ■ 2.8tt. For convenience all units 
have been scaled such that /i = 1. 

The parameter C was adjusted to reduce the errors due to the finite cutoff 
Since the spectrum of H is bounded below, errors 
generally drive the estimated eigenvalues higher. 
One strategy, therefore, is to optimize C by minimizing the trace of H restricted 
to the subspace Vk,-^^^,d„^i,^- The approach used here is a slight variation 
of this — we have minimized the trace of H restricted to a smaller subspace 
Vk' ,d' C Vk^^^,d^^^- The aim is to accelerate the convergence of the 
lowest energy states rather than thc entire space Vk^.^^.d^.^^- Throughout our 
analysis we used K'^^^, D'^^^ = 8,3. 

For |j = 0.50 we have plotted the imaginary part of A in Figure 1 and the 
real part of Amp in Figure 2. The value ^ = 0.50 is above the threshold for 
reliable perturbative approximation^ [^ < 0.25) but below the critical value at 
which (j) —> —(j) symmetry breaks spontaneously (^ w 2.5). The contribution of 
the one-particle state appears near p^ = (0.93)^ and the three-particle threshold 
is at approximately p^ — (2.9)^. We have chosen several different values for 
-^max, ^max, i'max to show the Convergence as these parameters become large. 
The plot for A^max, ATmax, -Dmax — 9, 19, 7 appears relatively close to the asymp- 
totic limit. The somewhat bumpy texture of the curves is due to the finite size 
of our periodic box and diminishes with increasing L. From dimensional power 
counting, we expect errors for finite A^max to scale as A^max- Assuming that 
-ft^max and -Dmax also function as uniform energy cutoffs, we expect a similar 
error dependence - and it appears plausible from the results in Figures 1 and 2. 
A more systematic analysis of the errors and extrapolation methods for finite 
A^max, ATmax, ^max, and L, wiU bc discussed in future work. 
>L 



In Figures 3 and 4 we have compared our spectral calculations with the 
two-loop perturbative result for ^ = 0.01. We have used iVmax, K^aax, -Dmax = 
9, 19, 7, and the agreement appears good. For small ^ the propagator has a 
very prominent logarithmic cusp at the three-particle threshold, which can be 
seen clearly in Figures 3 and 4. 

In Figures 5 and 6 we have compared results for |y ~ 0.25, 0.50, 1.00. We 
have again used iVmax, ^max, i'max — 9, 19, 7. In contrast with the quadratic 
scaling in the perturbative regime, the results here scale approximately linearly 
with jf . An interesting and perhaps related observation is that the magnitude 
of the multi-particle contribution to A remains small (< 0.003) even for the 
rather large coupling value |f = 1.00. 

4 Limitations and new ideas 

We now address the computational limits of the diagonalization techniques pre- 
sented in this work. These techniques are rather straightforward and can in 
principle be generalized to any field theory. In practise however the Fock space 
^ifmax-i'max bccomcs prohibitlvcly large, especially for higher dimensional theo- 
ries. The data in Figures 1 and 2 and crosschecks with Euclidean Monte Carlo 
results^ suggest that for A^max — 9 and L = 2.07r, • • • 2.87r our spectral results 
with ATniax, -Dniax = 19, 7 and A < i are within 20% of the Kmax, -Dmax -^ oo 
limit. In this case Vk^^^,d^^^ is a 2036 dimensional space and requires about 
100 MB of RAM using general (dense) matrix methods. 

Sparse matrix techniques such as the Lanczos or Arnoldi schemes allow us 
to push the dimension of the Fock space to about 10^ states. This may be 
sufficient to do accurate calculations near the critical point ^ ~ 2.5 for larger 
values of L and A^max- It is, however, near the upper limit of what is possible 
using current computer technology and existing algorithms. Unfortunately 
field theories in 2 -|- 1 and 3-1-1 dimensions will require much larger Fock spaces, 
probably at least 10"'^^ and 10^* states respectively. In order to tackle these 
larger Fock spaces it is necessary to venture beyond standard diagonalization 
approaches. The problem of large Fock spaces (^^10^ states) is beyond the 
intended scope of this analysis. But since it is of central importance to the 
diagonalization approach to field theory we would like to briefly comment on 
current work being done which may resolve many of the difficulties. The new 
approach takes advantage of the sparsity of the Fock-space Hamiltonian and the 
approximate (quasi-)sparsity of the eigenvectors. A detailed description will be 
provided in a future publication |9| . 

We start with some observations about the eigenvectors of the 4>f^i Hamil- 
tonian for A^max = 9, L = 2.57r and ATmax, -Dmax = 19, 7. To make certain that 
we are probing non-perturbative physics we will set ^ = 2.5, the approximate 
critical point value. We label the normalized eigenvectors as \vq), \vi),--- , 
ascending in order with respect to energy. We also define \bo), \bi),- ■ ■ as the 

^See pi for a discussion of these methods. 



normalized eigenvectors of the free, non-interacting theory. For any Vi we know 



Ek^^- 



1. 



(13) 



Let us define |l|'Ui)|L^i as the partial sum 



(") 



E \iKh)\'^ 



(14) 



fc=i. 



where the inner products have been sorted from largest to smallest 

\{b,,\v,)\>\{b,Jv,)\>--- . (15) 

Table 1 shows |||wi)||(„') for several eigenvectors and different values of n. 

Table 1 



\\\yi)\\M 


n = 10 


n = 20 


n = 40 


n = 80 


\vo} 


0.75 


0.84 


0.90 


0.94 


\yi) 


0.89 


0.92 


0.95 


0.97 


\V5) 


0.87 


0.91 


0.94 


0.96 


kio) 


0.77 


0.86 


0.90 


0.94 



Despite the non-perturbative coupling and complex phenomena associated with 
the phase transition, we see from the table that each of the eigenvectors can be 
approximated by just a small number of its largest Fock-space components. We 
recall that the Fock space for this system has 2036 dimensions. The eigenvec- 
tors are therefore quasi-sparse in this space, a consequence of the sparsity of 
the Hamiltonian. If we write the Hamiltonian as a matrix in the free Fock- 
space basis, a typical row or column contains only about 200 non-zero entries, 
a number we refer to as ./Vtransition- In H we show that a typical eigenvector 
will be dominated by the largest V -^transition elements.[| The key point is that 
\/-^transition is quitc manageable — on the order of 10'^ and 10^ for 2 + 1 and 
3-1-1 dimensional field theories respectively. Although the size of the Fock space 
for these systems are enormous, the extreme sparsity of the Hamiltonian sug- 
gests that the eigenvectors can be approximated using current computational 
resources. 

With this simplification, the task is to find the important basis states for 
a given eigenvector. Since the important basis states for one eigenvector are 
generally different from that of another, each eigenvector is determined indepen- 
dently. This provides a starting point for parallelization, and many eigenvectors 
can be determined at the same time using massively parallel computers. In 
we present a simple stochastic algorithm where the exact eigenvectors act as 
stable fixed points of the update process. 

^There are some special exceptions to this rule and they are discussed in [BJ. But these 
are typically not relevant for the lower energy eigenstates of a quantum field Hamiltonian. 



5 Summary 

We have introduced a spectral approach to periodic field theory and used it 
to calculate the propagator in 1 + 1 dimensional (f)"^ theory. We find that the 
straightforward application of these methods with existing computer technology 
can be useful for describing the multi-particle properties of the theory, infor- 
mation difficult to obtain using Euclidean Monte Carlo methods. However 
the extension to higher dimensional theories is made difficult by the large size 
of the corresponding Fock space. As a possible solution to this problem, we 
note that each eigenvector of the (1)1+1 Hamiltonian can be well-approximated 
using relatively few components and discuss some current work on quasi-sparse 
eigenvector methods. 
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Figures 

Figure 1 : Imaginary part oi A{p^) for ^ ~ 0.50 and several values for N^ax , -K^max , D^ 

Figure 2: Real part of Amp(p^) for ^ = 0.50 and several values for iVmax, ^max, -Dmax- 

Figure 3: Imaginary part of A(p^) for -rj = 0.01 and comparison with the two-loop 
result. 

Figure 4: Real part of Amp(p^) for |j — 0.01 and comparison with the two- loop 
result. 

Figure 5: Imaginary part of A(p^) for ^j- = 0.25, 0.50, 1.00. 
Figure 6: Real part of A^p{p^) for ^ = 0.25,0.50, 1.00. 
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Figure 3 
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Figure 5 



X/A\ 


= 0.25 


X/A\ 


= 0.50 


■--- X/A\ 


= 1.00 




P^ 



Figure 6 



Re[A^^(p2)] N„-,^,K„,^,D 



-Tnp 



0.003- 



0.002 



0.001 - 



max' -^^ax' '^max 



9,19,7 



-0.001 - 






X/A\ 


= 0.25 


— X/A\ 


= 0.50 


■--- X/A\ 


= 1.00 




P^ 



.'\.->-- 



12 



